Visualization of Coherent Destruction of Tunneling in an Optical Double Well System 
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We report on a direct visualization of coherent destruction of tunneling (CDT) of light waves in 
a double well system which provides an optical analog of quantum CDT as originally proposed by 
Grossmann, Dittrich, Jung, and Hanggi [Phys. Rev. Lett. 67, 516 (1991)]. The driven double 
well, realized by two periodically-curved waveguides in an Er:Yb-doped glass, is designed so that 
spatial light propagation exactly mimics the coherent space-time dynamics of matter waves in a 
driven double-well potential governed by the Schrodinger equation. The fluorescence of Er ions is 
exploited to image the spatial evolution of light in the two wells, clearly demonstrating suppression 
of light tunneling for special ratios between frequency and amplitude of the driving field. 

PACS numbers: 42. 50. Hz, 03.65.Xp, 42.82.Et 



Control of quantum tunneling by external driving fields 
is a subject of major relevance in different areas of physics 
p], Q. The driven double- well potential has provided 
since more than one decade a paradigmatic model to in- 
vestigate tunneling control in such diverse physical sys- 
tems as cold atoms in optical traps, superconducting 
quantum interference devices, multi-quantum dots and 
spin systems. Depending on the strength and frequency 
of the driving field, suppression 0, 0] or enhancement [H| 
of tunneling can be achieved. Tunneling enhancement 
is usually observed for high field strengths and driving 
frequencies close to the classical oscillation frequency at 
the bottom of each well. Since the enhancement gener- 
ally involves a transition through an intermediate state 
which is chaotic for strong enough driving amplitudes.it 
is often referred to as "chaos-assisted tunneling" p], Q. 
Observations of chaos-assisted tunneling have been re- 
ported in atom optics experiments [7] and in electromag- 
netic analogs of quantum mechanical tunneling [8|, uj • In 
particular, tunneling enhancement has been observed in 
two coupled optical waveguides [9[. In the opposite limit, 
Grossmann, Hanggi and coworkers [3[ found that, for cer- 
tain parameter ratios between amplitude and frequency 
of the driving, tunneling can be brought to a standstill. 
They termed this effect "coherent destruction of tunnel- 
ing" (CDT) and, since then, it has been of continuing 
interest. Driven tunneling is related to the problem of pe- 
riodic nonadiabatic level crossing and Landau-Zener (LZ) 
transitions p], [Io[. In particular, in the strong modula- 
tion limit CDT may be viewed as a destructive interfer- 
ence effect [Toj| . In spite of the great amount of theoretical 
work devoted to CDT, to date most of experimental evi- 
dences of CDT are rather indirect. In condensed- matter 
systems, dephasing and many-particle effects make tun- 
neling control more involved [11]. In Ref.[12] coherent 
control of Rabi oscillations in Josephson-j unction circuits 
irradiated by microwaves has been reported, however the 
condition for CDT was not reached. Quantum inter- 
ference effects and evidences of CDT in qubit systems 
have been recently reported in pjj QjJ , whereas suppres- 
sion of quantum diffusion, also known as dynamic local- 



ization, has been experimentally demonstrated in Refs. 
PH EH, E3 • However, CDT is a rather distinct effect than 
dynamic localization (see [H.[l8l|). For a cleancut demon- 
stration of CDT, a direct visualization of the dynamics 
is desirable, which was not accomplished in all these pre- 
vious experiments. Engineered optical structures, on the 
other hand, have been recently demonstrated to provide 
a very appealing laboratory tool for a direct visualiza- 
tion of optical analogs of quantum mechanical phenom- 
ena which require a high degree of coherence [20[ . 
In this Letter we report on the first visualization of 
CDT dynamics using an optical analog of a driven 
bistable Hamiltonian based on two tunneling-coupled 
curved waveguides [19[ which enables an experimental 
access to the full space-time evolution of the corre- 
sponding quantum mechanical problem Q. The struc- 
ture designed to visualize CDT consists of a set of two 
L = 24-mm-long parallel channel waveguides, placed at 
a distance a = 11 //m, whose axis is sinusoidally bent 
along the propagation distance z with a bending pro- 
file xq(z) = Acos(2ttz/ A) [see Fig. 1(a)]. The waveguides 
have been manufactured by the ion-exchange technique 
(2lj in an active Er-Yb phosphate substrate and probed 
at A ^ 980 nm wavelength using a fiber-coupled semi- 
conductor laser [Fig. 1(b)] with ~ 8 fim mode diameter. 
A transverse scan of the fiber along the sample is used 
to preferentially excite either one of the two wells. The 
probing light is partially absorbed by the Yb 3+ ions (ab- 
sorption length ~ 6 mm), yielding a green upconversion 
luminescence arising from the radiative decay of higher- 
lying energy levels of Er 3+ ions (22[. By recording, at 
successive propagation lengths, the fluorescence from the 
top of the sample using a CCD camera connected to a mi- 
croscope (magnification factor ~ 12) mounted on a PC- 
controlled micropositioning system, we could trace with 
accuracy the flow of light along the sample. It was previ- 
ously shown [19[ that the evolution of light waves in the 
optical double-well system is formally equivalent to the 
dynamics of a periodically-driven nonrelativistic quan- 
tum particle in a double-well potential. In the scalar and 
paraxial approximations, light propagation is described 
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FIG. 1: Microscope images (top view) of the sample showing 

(a) a few set of manufactured coupled optical waveguides, and 

(b) the fiber coupling geometry for waveguide excitation. 



by the equation @, [l9[ 



reduced wave- 
~ n s - n(x,y) 



where: A = A/(27r) = 1/k is the 
length, V"0,2/) = [n 2 s - n 2 (x, y)]/(2n s ) 
is the double-well potential, n(x,y) is the refrac- 
tive index profile of the two-waveguide system, 
and n s is the reference (substrate) refractive index. 
The quantum-optical analogy can be retrieved af- 
ter a Kramers-Henneberger transformation [l?], EH 
x' = x - x (z), y' = y, z' = z, <t>{x' \y' ', z') = 

tj){x f , y', z') exp -i(n s /X)x (z')x f - i(n s /2X) f* d£ ) 

where the dot indicates the derivative with respect to z' \ 
and after elimination of the ^/-dependence of the field (j) 
using a standard effective index method [23[. Equation 
(1) is then transformed into the following Schrodinger 
equation for a particle of mass n s in the double-well 
potential V e {x r ) ~ n s — n e (x r ) under the action of a 
sinusoidal force F(z') [19j : 



iX 



dz' 



X^d^ 
2n s dx' 



r+Veix'W-Fx't^Tiot-Fx'fa (2) 



where n e (x') is the effective index profile of 
the waveguide system and F(z f ) = —n s xo(z f ) — 
(47r 2 An s /A 2 ) cos(27rz'/A) is the ac force. Note that, in 
the optical analog, the Planck constant is played by the 
reduced wavelength A of photons, whereas the tempo- 
ral variable of the quantum problem is mapped into the 
spatial propagation coordinate z' . CDT is thus simply 
observed as a suppression of photon tunneling between 
the two waveguides along the propagation direction. 
The refractive index profile n(x,y) has been measured 
by a refracted-near-field profilometer (Rinck Elektronik) 
at 670 nm. The measured 2D index profile is depicted 
in Fig. 2 (a), together with its section profiles along the 
horizontal (HH) and vertical (VV) lines. Figure 2(b) 
shows the corresponding symmetric double-well poten- 
tial V e (x) ~ n s — n e (x) obtained by the effective index 
approximation. To derive n e (x), the measured 2D index 
profile n(x,y) was fitted [see dashed curves in Fig. 2(a)] 
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FIG. 2: (color online) (a) Measured refractive index profile of 
the double waveguide system, (b) Computed profile of the ef- 
fective ID double-well potential, (c) Manifold of quasi-energy 
crossing in the (A, A) plane (solid line). The dashed line has 
been evaluated from the first zero of the Bessel function using 
Eq.(3). The inset is an enlargement of the linear portion of 
the manifold, corresponding to CDT. The five points in the 
figure correspond to the geometrical parameters of the optical 
waveguides manufactured in our experiment. 



by the relation [24[ n(x, y) ~ n s + An[g(x — a/2) + g(x + 
a/2))f(y), where An c± 0.0124 is peak index change, 
g(x) = [eYf((x + w)/D x )-eif((x-w)/D x )}/[2eif(w/D x )} 
and f(y) = exp(—y/D y ) define the shape of the in- 
dex profile parallel to the surface of the waveguide (in- 
direction) and perpendicular to the surface (^-direction), 
respectively, 2w c± 5 jam is the channel width and D x ~ 
4.3 /im, D y c± 3.3 fim are the lateral and in-depth dif- 
fusion lengths. The numerical computation of the eigen- 
values for the Hamiltonian Ho = —X 2 /(2n s )d 2 , + V e (x f ) 
in absence of the driving force indicates that the double- 
well potential V e supports four bound modes £,i(x') whose 
energies F\ (I = 1,2,3,4) are depicted in Fig. 2(b) by 
the four horizontal solid lines. The linear combinations 
ur,l{x') = [£,i(x') ± ^2(x')]/V2 of the eigenfunctions 
£,i(x') and £2(2/) associated to the quasi-degenerate en- 
ergy levels Ei and E2 below the barrier correspond to 
photon localization in the right (R) or in left (L) well 
of the potential, so that an initial excitation of one of 
the two wells, obtained by launching the light into either 
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FIG. 3: Fluorescence light distribution (top view) in the 
straight waveguide coupler, and corresponding light intensity 
distribution predicted by Eq.(l), for different coupling condi- 
tions, (a) Excitation of one of the two waveguides; (b) Beam 
launching in the middle between the two waveguides. 



one of the two waveguides, is given approximately by the 
superposition of the two lowest eigenstates £i(V) and 
<^2(V)- For the straight waveguides, the field evolution 
is dominated by the splitting of this doublet, leading to 
a periodic tunneling of photons between the two waveg- 
uides with a spatial period d\2 = 2ttX/(E2 — E\) ~ 7.94 
mm. This is clearly shown in Fig. 3 (a), where the mea- 
sured fluorescence pattern corresponding to initial ex- 
citation of one of the two waveguides is reported. For 
the sake of clearness, in the picture the luminosity level 
of the fluorescence, which decreases with propagation 
distance due to light absorption, has been gradually 
rescaled at successive frames. The measured pattern 
is very well reproduced by a direct numerical simula- 
tion of Eq.(l), performed with a standard beam prop- 
agation software (BeamPROP 4.0). Higher-order bound 
modes of the double- well potential shown in Fig. 2(b) 
can be excited by different beam launching conditions. 
For instance, if the fiber is positioned in the middle be- 
tween the two waveguides, an excitation of the modes 
£i (x) and £3 (x) with even symmetry is preferentially at- 
tained. In this case, the field evolution is dominated 
by the splitting of the states £i(x) and £3(2;), leading 
to a periodic fluorescence pattern with the short spatial 
period d 13 = 2ttX/(E 3 - Ex) ~ 545 /im [see Fig.3(b)]. 
The tunneling dynamics in presence of the external force 
F strongly depends on the amplitude and (spatial) fre- 
quency uo = 27r/A of the force as compared to the en- 
ergy level spacing of the double- well system p|. For in- 
stance, at modulation frequencies comparable with the 
frequency spacing (£3 — E2) /A, tunneling is expected to 
be enhanced. This case was previously demonstrated for 
a two- waveguide optical system in Ref.[9j]. Conversely, 
CDT occurs approximately for a modulation frequency 
in the range (E2 — E\) < Xuo < (E3 — E2) and for a mod- 
ulation amplitude which corresponds to exact crossing 
between the quasienergies e\ and 62 associated to the low- 
est tunnel doublet [1]. In Fig. 2(c) the solid curve shows 
the manifold associated to the exact crossing €2 = e\ in 
the (A, A) plane, as numerically computed by means of 
a two-level approximation of the related driven tunnel- 
ing problem Q. In the high- frequency limit, i.e. for 
Xuo ^> (E2 — Ei) but Xuj < (E3 — E2) to avoid the partic- 
ipation in the dynamics of the third level of energy £3, 
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FIG. 4: Measured fluorescence light distribution in the 
"waveguide" reference frame (x,z) for the straight waveguide 
coupler (a), and for the three curved waveguide couplers [(b), 
(c) and (d)] with increasing period A and amplitude A corre- 
sponding to points 1,2 and 3 of Fig. 2(c). 
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FIG. 5: Measured fluorescence light pattern, and corre- 
sponding photon density \tp(x, 0, z) | 2 predicted by Eq.(l), in 
curved waveguide couplers corresponding to (a) break-up of 
quasienergy crossing [point 4 in Fig. 2(c)], and (b) stroboscopic 
suppression of tunneling [point 5 in Fig. 2(c)]. 

an approximate expression for the quasienergy difference 
Ae = e 2 -ei reads 0,11, EH 

Ae=(E 2 -E 1 )J {^^) (3) 

where \i\2 = (£i|#|£2)- Considering the first zero of the 
Bessel function, the condition Ae = is thus represented 
by a straight line in the (A, A) plane, which is depicted by 
the dashed curve in Fig. 2(c). Such a curve, however, de- 
viates from the exact one as A increases and approaches 
47rA/(£^2 — Fi) ~ 16 mm, where the solid curve drops 
toward zero. Crossing of the quasienergies is a necessary 
-but not a sufficient- condition for the occurrence of CDT. 
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In fact, CDT requires additionally that the degenerate 
Floquet states at exact energy crossing do not show ap- 
preciable amplitude oscillations in one period. A detailed 
numerical analysis of Eq.(l) shows that CDT indeed oc- 
curs along the linear portion of the manifold of Fig. 2(c), 
i.e. for A < 2.5 mm, which is represented by the enlarged 
inset in the figure. We experimentally demonstrated the 
onset of CDT in this portion of the manifold by manu- 
facturing three curved waveguide couplers corresponding 
to the points 1,2 and 3 of Fig. 2(c), i.e. to A = 1.5 mm, 
A ~ 5.5 jam (point 1), A = 2 mm, A ~ 7.3 jam (point 2) 
and A = 2.4 mm, A ~ 8.8 jam (point 3). Figure 4 shows 
the measured fluorescence patterns, as recorded on the 
CCD camera from the top of the sample, for the straight 
waveguide coupler and for the three curved waveguide 
couplers when one of the two waveguides is excited at 
the input plane. Note that, since the fluorescence is 
proportional to the local photon density, the patterns 
in the figure map the profile of |^| 2 in the "waveguide" 
reference frame. Therefore, in Figs. 4(b), (c) and (d) 
the condition for CDT is clearly demonstrated because 
the photon density follows the sinusoidal bending profile 
xo(z) of the initially excited well, without tunneling into 
the adjacent waveguide. We also experimentally checked 
that the observation of CDT requires that the follow- 
ing two conditions must be simultaneously satisfied: (i) 
quasienergy crossing, and (ii) absence of appreciable am- 
plitude oscillations of the degenerate Floquet eigenstates 
within one oscillation period [2, H|. As an example, in 
Fig. 5 (a) we show the measured fluorescence pattern -and 



corresponding photon density pattern predicted by the 
theory- for the curved waveguide coupler corresponding 
to point 4 of Fig. 2(c) (A = 2.5 mm and A — 7.9 jam), 
in which the condition (i) is not satisfied. Note that 
in this case the pattern periodicity is broken and tun- 
neling is not suppressed, though the tunneling rate is 
reduced as compared to the straight waveguide coupler 
case [compare Fig. 4(a) and Fig. 5(a)]. Figure 5(b) shows 
the measured fluorescence pattern for the curved waveg- 
uide coupler corresponding to point 5 of Fig. 2(c) (A = 6 
mm and A ~ 20.5 /am). In this case the condition (i) for 
quasienergy crossing is fulfilled, however over one oscilla- 
tion period the Floquet eigenstates show non-negligible 
amplitude oscillations. Though tunneling is suppressed 
at the stroboscopic distances z = A, 2 A, over one os- 
cillation period an appreciable fraction of light is ob- 
served to tunnel forth and back between the two waveg- 
uides. Such a stroboscopic destruction of tunneling can 
be viewed as a result of destructive interference between 
successive LZ transitions taking place at periodic level 
crossings p], [Io|, i.e. at the positions z = A/4,3A/4, ... 
where xq = 0. This periodic regime, however, does not 
correspond to a true CDT 

In conclusion, we reported on the first visualization of 
photonic CDT in an optical double-well system which 
mimics the corresponding quantum-mechanical problem 
originally proposed in [3[. The two basic conditions for 
the observation of CDT, namely quasienergy crossing and 
absence of amplitude oscillations of the degenerate Flo- 
quet doublet, have been experimentally demonstrated. 
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